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We introduce a deterministic model defined on a two dimensional hyper- 
bolic lattice. This model provides an example of a non random system whose 
multifractal behaviour has a number theoretic origin. We determine the mul- 
tifractal exponents, discuss the termination of multifractality and conjecture 
the geometric origin of the multifractal behavior in Liouville quasi-classical 
field theory. 



I. INTRODUCTION 

The concept of multifractality consists in a scale dependence of critical exponents [|l|. 
It has been widely discussed in the literature in the context of various problems such as, 



for example, statistics of strange sets 0-01 13 1, diffusion limited aggregation M, wavelet 



transforms |^ , conformal invariance . This concept also proves to be useful in the context 



of disordered systems [^,^. It was recently found that the ground state wave function 
of two dimensional Dirac fermions in a random magnetic field has a multifractal behavior. 
The field theoretic investigation of the multifractality has been undertaken in the papers , 
while different interpretations of these field theoretic results from a geometrical and physical 
points of view were presented in B and [IlOl correspondingly. This problem was recently 



reanalyzed in the more general setting of systems caracterized by logarithmic correlations 



Our work is mainly inspired by the approach developed in where the authors obtain 
the multifractal exponents of the critical wave function by a mapping on the problem of 
directed polymers on a Cayley tree. However our starting point is different and we treat a 
deterministic model defined on a Cayley tree. We take advantage of the fact that the Cayley 
tree can be isometrically embedded in a space of constant negative curvature. We assume 
that each vertex of the tree carries a Boltzmann weight that depends on the hyperbolic 
distance from a given root point. The corresponding partition function is a sum over a 
finite number of tree vertices and has the form of a truncated Poincare series. Its scaling 
dependence on the size of the system is controlled by the probability distribution of traces 
of 2 X 2 matrices which belong to a discrete subgroup of P5'L(2,IR). This distribution. 



obtained by using the central limit theorem for Markov multiplicative processes |30], allows 
us to compute the multifractal exponents and discuss the termination of multifractality. 
The study of the convergence of the measure on the boundary reveals some interesting links 



with work of Gutzwiller and Mandelbrot 13] on multifractal measures. Another interesting, 
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although more speculative, aspect is connected with a geometric approach to Liouville field 
theory arising in the study of low dimensional disordered systems [3,Pl p8i^ . We suggest 
that in two dimensions our model exhibits a new type of multifractal behavior which has a 
purely geometric origin. 

This paper is organized as follows. In section we introduce the geometrical model 
possessing the multifractal behavior, develop methods for its investigation and explicitly 
show the number theoretic origin of multifractality; section |Tl| is devoted to applications 
of these results to quasi-classical 2D Liouville field theory (LFT); the conclusion presents 
some speculations regarding the applicability of our geometric considerations to some other 
disordered physical systems. 



II. THE MODEL 



We begin with the investigation of geometrical properties of lattices uniformly embedded 
in the hyperbolic 2-space. Lattices under consideration are defined as follows: we construct 
the set of all possible orbits of a given root point under the action of a discrete subgroup 
of PSL{2, IR) (group of motion of the hyperbolic 2-space). We restrict ourselves with the 
simplest example of 3-branching Bethe lattice (Cayley tree) which is generated by refiections 
of zero-angled curvilinear triangle — see fig.0. 




FIG. 1. A Cayley tree in the Poincare disc. Sample points: A(-^j, {5, C} ( =F^ - y 




The graph connecting the centers of the neighboring triangles forms a Cayley tree iso- 
metrically embedded in the Poincare unit disc (a Riemann surface of constant negative 
curvature) . 

Consider the n^^ generation of the vertices of the 3-branching Cayley tree. Denote by 
rj{n) the Euclidean distance of the vertex j (which belongs to the n^^ generation of the 
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Cayley tree) from the center of the unit disc (1 < J < 3 x 2" ^). The corresponding 
hyperbohc (geodesic) distance dj{n) is given by: 

"'("J^'-T^ (.,(«) = tani^) (1) 

Define the generating function Z{q, N) 

N /3x2"-i \ 

Z{q,N)^j:{ E (2) 

n=l \ j=l / 

In a physical context Z{q, N) may be interpreted as a partition function on the hyperbohc 
lattice with an action linear in the length of trajectory. 

The Bethe lattice involved can be constructed by the action of the discrete group Te 
which operates on the unit disc by a set of fractional-linear transformations. Despite the 
simple structure of the group it is believed that the techniques involved are quite general 
and could be easily generalized in order to cover more sophisticated lattices. 

We are interested in the scaling dependence of the partition function Z{q, N) as a function 
of the size N of the system. Scaling considerations suggest the following behaviour 

Z{q, N) = (3) 

where in our case £max = 3(2^ — 1) is the total number of Cayley tree vertices in the bulk 
restricted by the generation n = N. 

We show below that the critical exponent T{q) defined as follows 

lnZ{q,N) . , 

depends nonlinearly on q i.e. exhibits the multicritical behavior (note that r* = r/ln2). 
Note that the free energy normalized per volume of the system f(q, N) = — ''^-^fa-'^) coincides 
with the multifractal exponent r(g): 

hm f{q,N)^T{q) (5) 

AT— »oo 



A. Numerical results 



We first compute numerically the histogram, which counts the number of vertices be- 
longing to generation n (properly normalized), Wn{d), lying in the shell [d, d + Sd]. 

In our particular computations we restrict ourselves with two cases depending on the 
length of the trajectories: 



3 



1. Short trajectories. We enumerate all trajectories and the computations have been 
carried out for all n e [1, iV] up to = 25 generations. The figure fig.^a shows the 
histogram for the distribution of hyperbolic distances for n = 25. The absolute value 
of number of events in the fig.@a depends on the particular choice of the width of the 
shell 6. It can be seen from fig.^a that the corresponding plot is highly nonsymmetric 
with respect to the mean value (d). 

2. Long trajectories. For n = 200 the enumeration of all different paths is very time 
consuming, therefore we compute numerically the histogram Wn{d) developing partial 
ensemble of 200 000 directed random walks of n = 200 step each. As n ^ cxd the 
distribution function Wn{d) becomes more and more symmetric in accordance with 
the statement that there exists a central limit theorem for such random walks on 
noncommutative groups (see the discussion below). The results of corresponding nu- 
merical computations are presented in Fig.||b. The distribution Wn{d) is well fitted by 
a Gaussian function: 

(d-<d))2 

Wn{d) = Aoe-^^ 

where for n = 200 one has: Ao ^ 1929.96 and depends on normalization; (d) ~ 159.18; 
A2 ^ 17.01. 




FIG. 2. Distributions Wn{d) up to normalization, compared to their Gaussian fits: one can 
notice the slow convergence from strongly nonsymetric regime for n = 25 (a) to a Gaussian regime 
for n = 200 predicted by the central limit theorem (b). 



In spite of the fact that convergence to the Gaussian distribution is slow, the linear 
dependence in n of the mean value (d) = 7n and the variance {{d — (rf))^) = = a^n is 
numerically evident, which permits one to get an accurate estimate of 7 and a^. 
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The numerical computation of the probabihty distribution Wn{d) allows one to compute 
the multifractal exponent r(g) following the definitions The corresponding results 

are shown in fig.^, for iV = 40. Due to the slow convergence of the distribution, the 
discrepancy between numerical data (technically limited to < 40) and the theoretical 
prediction can not be quantitatively taken into account. We here insist on the multifractal 
behaviour, shown by the non-linear depence on q. 
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FIG. 3. Multicritical behavior of the exponent T{q) for N = 40, compared to theoretical predic- 
tion. 



B. Analytic results 



Let us return to the definition of the model and recall that the group Tq acts in the 
hyperbolic Poincare upper half-plane Ti.^ = {z E C, lm(2;) > 0} by fractional-linear trans- 
formsQ. The matrix representation of the generators of the group Tq is well known (see. 



for example |]23|), however for our purposes it is more convenient to take a framework that 
consists of the composition of standard fractional-linear transform and complex conjugacy. 
Namely, denoting by z the complex conjugate of z, we consider the following action 

az + b 

(6) 




cz + d 



A possible set of generators is then: 



It is convenient first to define the representation of the group Tg in the Poincare upper half-plane 
and then use the conformal transform to the unit disc. 
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1 -2/v^\ . _(l 2/v^\ 
-1 j' 1^0 -1 j 



1/^3' 
73 



(7) 



Choosing the point (a;o, iyo) = (0, i) as the tree root — see fig.§, any vertex on the tree is 

n 

associated with an element M„ = JJ^ where G {0, 1, 2} and is parametrized by its 

fc=i 

complex coordinates Zn = M„((— in the hyperbolic plane. 




FIG. 4. Poincare hyperbolic upper half-plane H?' 

Strictly speaking should be identified with S'L(2, IR)/S'0(2); we here identify an 
element with its class of equivalence of S0{2). If one denotes by d{Mn) = d{i,Zn) the 
hyperbolic distance between i and Zn, the following identity holds 



2 cosh (rf(M„)) = Tr(M„Mt) 



(8) 



where dagger denotes transposition. 



1. Distribution function, invariant measure on the boundary and Lyapunov exponents 



We are interested in the distribution function Wn{d) which is the probability to find the 
tree vertices in generation n at the distance d from the root point. It means that we are 
looking for the distribution of the traces for matrices M„ which are the irreducible products 
of n generators. If we denote by l{Mn) the irreducible length of the word represented by the 
matrix M„, then Mn is irreducible if and only if l{Mn) = n. Such word enumeration problem 
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is simple in case of the group Tg, because of its free product structure: Fe ~ Z2 ® Z2 ® Z2. 

n 

Indeed, if M„ = JJ^ h^^ one has l{Mn) = n if and only if h^^ 7^ ^a-k-i Hence we have to 

k=l 

study the behavior of the random matrix M„, generated by the following Markovian process 

{(a„ + 1) mods with probability | 
[On + 2) mods with probabihty ^ 

We use the standard methods of random matrices and consider the entries of the 2x2- 
matrix M„ as a 4- vector V„. The transformation M„_|_i = M„/iq, reads 



^n+l = .t (10) 




This block-diagonal form allows one to restrict ourselves to the study of one of two 2- 
vectors, composing V„, say f„. Parametrizing t>„ = (^)„ cos6'„, sin 6'„) and using the re- 
lation d{Mn) = dn — 2\ia Qn Valid for n 1, one gets a recursion relation Vn+i = hl^Vn in 
terms of hyperbolic distance dn- 



dn+i = dn + ln 



5 4 

- + - cos(26'„ + v^o 



(11) 



where (pa depends on the transform ha through = (2a — l)vr/3 (a = 0, 1,2), while for 
the angles one gets straightforwardly 

tan6'„+i = /ia(tan(6'„)) (12) 

and the change a — > a means the substitution (0,1,2) — » (1,0,2). Action of ha is still 
fractional-linear. 

Define now three invariant measures /ia(^) corresponding to transformations of M„ (n 3> 
1) whose last step is given by a matrix ha- The form of (|l^) suggests to consider the 
corresponding fia{x) with x = ta.n9. We are then led to study the action of Tg restricted on 
the real line parametrized by x. Interesting properties of the average JI = (/io + /ii + /i2)/3 
have been discussed by Gutzwiller and Mandelbrot [|T^. In particular they pointed out the 
connexion with the arithmetic function /3(^) which maps some number ^ G [0, 1] written as 
a continued fraction expansion 

1 

1 

ni + 



n2 + . . . 

to the real number (3 whose binary expansion is made by the sequence of ni — 1 times 0, 
followed by 77-2 times 1, then times 0, and so on. To account for this fact, one has first 
to notice that the construction (^) of any word M in Tq is exactly encoded by the binary 
representation of a real ^, the n's letter of this expansion being — a„ + ImodS. The 
second argument, due to Series [Q, is that the real part of the vertex M{i) is precisely the 
continued fraction ^{C)- Therefore (3{^) has to be proportional to the "number" of vertices 
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lying in the interval [0,^], that is to /i([0,^]). Taking the limit n — > oo is not well defined. 
An alternative, which was used in this work, is to define the limit of the following 

recurrency: 

E (13) 



dha(x) 



dx 



The symetry of such expression leads, after summing over a, to the following relation ad- 
mitting as fixed point 7l(x) at n ^ cxd: 



Q=0 



dhn(x) 



dx 



(14) 




FIG. 5. Invariant measure /x as a function of 9 



The convergence fi 



for n — > oo is assured by ergodic properties of such 



functional transform in case of equation (14), and has been successfully checked numerically 
by comparing to direct sampling of different orbits. Obtaining fia for a = {0, 1, 2} from ]I 
is not difficult, taking into account the symmetric role that they play with respect to the 
three intervals Jq =] - oo,-l/V3], h = [-1/^3,1/^/3], h = [1/^3, +oo[ (see fig. prBII ). 
Contracting properties of ha{x) allow convergence of (O) only if 



/ifc(x) = 3x4(a;)yu(x) 



(15) 



where Xh characteristic function of the interval J^. 
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We would like to point out an interesting fact, even if far from being rigorous, which is 
very similar to the argument put forward in for justifying the connexion between the 
invariant measure and the arithmetic function It has been shown in |^2[ that the lattice 
under consideration can be isometrically embedded in a 2-manifold A4 = {c\ri{z)\^, z G Ti.^}, 
where 



ri{z) 



^TTiz/12 



na 



^2ninz\ 



n=l 



is the Dedekind //-function. The mountain range (relief) displays a very steep valley 
structure, and our tree lattice was defined as the ridges of this relief. The natural "counting" 



of vertices whose real part lies in [0, ^] in [13| is in our case equivalent to counting the number 
of maxima of \ri{x + iO"*")!^, that can be directly reexpressed as a density if one admits that 
all maxima are equivalent and well separated: 



\7]{x + i0+)\^ 



(16) 



The intriguing fact is that r]^ is an automorphic form of weight 2, what makes \ri\'^ precisely 
a possible fixed point of Eq.(0). We recall the fundamental property of automorphic forms 
/ of weight 2 under the action of SL{2, IR): 



/(^) 



^i4>{a,b,c,d) ^ /az + b' 
{cz + d)'^ \cz + d^ 



(17) 



The main problem is that the boundary behavior of automorphic forms is far from trivial 
(see [^), and (|16D has no rigorous mathematical sense. In particular compatibility of (p!4[) 
and ([T^) is not obvious even numerically. Nevertheless we insist on the fact that fik is defined 
with no ambiguity by ([T^), what enables us to compute the desired Wn{d). The crucial point 
here, already required for convergence of /i*^""-*, is ergodicity property of 9n- It means that 
for n 3> 1, the distribution of 9n is exactly given by 71(6'), independently of n and initial 
condition. Then, denoting by c?" the value dn obtained for a word ending with ha, one can 
transform ([TT|) in the following way: 



d 



n+l 



1 ^ 



rf" + ln 



5 4 

- + - cos(26'„ + ipi) 



(18) 



with the condition / ^ a. Thus we obtain 



1 ^ 

^ A:=0 jVfc 



/■V2 /5 4 

/ de^^ie) - + -cos(2^^ + y.,) 

J — 7r/2 \ O O 



5 4 



ik' 



(19) 



which finally leads to 



"^nid) = ^ r dke 

ZTX J-oo 



ikd 



d9fi,i9-^) (^^ + lcos2^ 



ik' 



(20) 
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This form suggests that for n large Wn{d) satisfies a central limit theorem. Indeed such 
a theorem exists (see [24,^|) for Markovian processes, provided that the phase space is 



ergodic. We are then led to compute only the first two moments (Lyapunov exponents) 
which gives 

7 = — = r d9^ii{e-^)\n(l + Uos2e) ^ 0.792 (21) 
n Jo b \6 6 J 

and 



a 



n 



12 -Y (22) 



with 

72 = Ij de^l^i9 - ^)\n' (^^ + ^cos29^ ^ 0.68 (23) 

Numerical simulations presented in previous section yield 7 ~ 0.793 and 72 ~ 0.66, which 
finally allow us to conclude that for n 1 Wn{d) has a Gaussian behavior 

Wn{d)=Ae-^^ (24) 

centered at 771 and of variance a^n [A is the normalization). 

The numerical values of the Lyapunov exponents 7 and 72 (see Eqs.(^) and (p3|)) are 
obtained by means of semi-numerical procedure which involves the numerical information 
about the invariant measure ^i{9). However one can get the estimates for the Lyapunov 
exponents 7 and 72 by approximating the measure [ii{9) on the interval < 6' < | in two 
different ways: 

^^,{9-l)^^,^{9) = - 

l 3 (2^) 

^^,{9--)^^^^{9) = -smm 

Both measures yU^ and yuf are properly normalized on the interval [0, |]. Substituting (^5[) 
in (|2T| ) and ( p3D and computing (analytically for 7) the corresponding integrals, one finally 
gets: 

0.665 ^1 • { 72 ^ 0.684 ^'^^^ 

As one can see, the agreement between numerical values of Lyapunov exponents obtained 
for the measures /^i and its approximants /if ''^ is reasonable for /if and very good for /if. 
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2. Multifractal exponents 



The partition function Z[q, N) introduced in (Q) can be defined for any discrete subgroup 
of PSL{2, IR) of generic element r by 

Z{q,N)= Yl e-^'^M (27) 

t,1{t)<N 

and the associated critical exponent is then 

The probabihty distribution ( P^ enables us to rewrite ( P7D for the group Vq in the limit 
n ^ 1 as follows 

N 

Z(g,iV) = ^3x2"^V (29) 

n=l 

where 

a„= / e-'i^Wn{t)dt= \ Ae'^^e-'i'dt (30) 
Jo Jo 

The following two cases should be distinguished: 

• For q < 7/0"^, the minimum of the exponent in Eg . (1301) is within the range of integra- 
tion and 

an ~ e-"^''+"'^''?'/2 (31) 

hence 

q N 

Z{q, iV) ~ - £ e'^(i"^-^'?+'^'^'/2) (32) 
2 fc=i 

The convergence of the sum (p2D for N 00 depends on the sign of the function in 
the exponent. For 

7- V7^-2a^ln2 

q<qo = (33) 

one has — In 2 + 75 — a'^q'^/2 > and the multifractal exponent is 

r(g) = -ln2 + 7g-aV/2 (34) 

while for q > qo, the series Z{q, N) converges and r(g) = what signals the termination 
of the multifractality. Note that go is real at least in the case of the group Tg. 
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For q > 7/0"^ the minimum of the exponent in Eq. (pOD is out of the range of integration 
and 

Inon--^ (35) 



hence Z{q, N) is no longer extensive in N, which leads to r(g) = r(go) = 0. 



The nonlinear dependence on q obtained above shows the multifractal behaviour of this 
model below the termination point go- It seems more transparent to summarize all these 
results in a table 



q < 






q > -f/a^ 


r{q) = r{qo) = 



3. Conformal mapping approach to computation of partition function and multifractal exponent 



We propose in this section a completely different approach allowing to get a closed 
analytic expression for the partition function similar to ^(g, A^'max) (see Eq.(0)). The con- 
struction presented below is a by-product of our former investigations of analytic structure 
of the covering Riemann space of the multi-punctured complex plane (see, for review [pT[]). 
We explore the properties of the Jacobian of conformal mapping of the infinite complex plane 
with a triangular lattice of punctures into the unit disc parametrized hj w = re 
this particular case represents the multi-sheeted universal covering space pl[ 
define two functions f{r,a) and g{r): 



which in 
Namely, we 



9[ (0, e''^' 



\l + iw\^ 



8/3 



= c 



\viaw))f 

\l + iw\^ 



9[r) 



1 



(36) 



where 



e[ (0,e*"^) =2e*^^^(-l)"(2n + l)e*""(" 



+i)C 



CM 



-in/3 



n=0 

w + e^'^l^ 



w — I 

-8/3 



(37) 



0.933293 



One can show that the functional equation 



-1 = (38) 

9\r) 
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has a family of solutions {tc, etc) exactly at positions of 3-branching Cayley tree isometrically 
embedded in the hyperbolic unit disc (in the Klein's model of the surface of constant negative 
curvature). In fig.^we have plotted the 3D section of the function 



u{r, a) 



f{r,a) 
9{r) 



(39) 



in polar coordinates (r, a) for 0.9 < u{r, a) < 1. The function u{r, a) has local maxima with 
one and the same value u = 1 only at the coordinates of isometric embedding of 3-branching 
Cayley tree in the hyperbolic unit disc. The proof of this fact is given in Appendix 0. 

Thus, we can rewrite (|^) in the following closed form (recall that Euclidean distance r 
and the corresponding hyperbolic distance d are linked by the relation (^) 



r{d) 2tt 



Z{q,d) = ^ J J e-^'"T^5(lnM(r,«)) 







d\nu{r, a) 



dr 



rdrda 



(40) 



where for (5[ln'u(x)] we use the standard integral representation (5[ln-u(x)] = ^ I b^i^)] 




FIG. 6. 3D parametric plot of the function u{r,a) in the section 0.9 < u{r,a) < 1. 



It is noteworthy to pay attention to the difference between the partition functions Z{q, N) 
(Eq.(]^)) and Z{q, d) (Eqs.(^) and (^DD). The function Z(g, A^) counts the weighted number 
of Cayley tree vertices up to the generation N for nonfixed maximal radius r{d) = tanh((i/2) 
in the hyperbolic unit disc, while the function Z{q, d) counts the weighted number of Cayley 
tree vertices within the hyperbolic disc of radius r{d) for nonfixed maximal generation A^. 
The last partition function is in fact related to the number of tree vertices inside the disc of 
radius d. This is the content of the famous circle problem first formulated by Gauss for the 
Euclidean lattice Z^. The extension to the non- Euclidean case is due to Delsarte |3l[] (see 
also 
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III. MULTIFRACTALITY IN 2D QUASI CLASSICAL LIOUVILLE FIELD 

THEORY 



Our starting point is the family of normalized wave functions ipki'^) defined as follows 

I Ifc -(/3(x) 

M^) = ^ T7? (41) 

where integration extends to a disc of radius R and the potential v'(x) is distributed with 
Gaussian distribution function 

P[¥.(x)] cx exp J cix(V^(x))^| (42) 

The problem defined in (^TD -(^) appears in various models which will be discussed in the 
next section. 

The multifractal exponent for the quenched and annealed distributions of disorder in 
(p])-p^) can be computed in the standard way 



Tquiq) = — lim - — — ^— for quenched disorder 

r^M = -Um'°"> for anuealed disorder 

R-*oo \nR 

where Q{q,R) = J (ix|-?/'fc(x)p'' and the brackets (...) denote averaging with the distribution 
^b(x)]. 

We pay attention to the case of annealed disorder and our aim is to evaluate the corre- 
lation function 

{Q{q,R))= |rfx(|^,(x)p^) (44) 

The averaging (...) in (^) means 

{■■■)sM = J ^^b(x)]e-^[-Wl (45) 

where 

^b(x)] = ^/rfx{(V^(x)f} (46) 

In order to take into account proper normalization of the wave function ipk{^) it is 
convenient to use a Lagrange multiplier A, so that eventually 

(Q(g, R)) = Jd^oJ P[(/;(x)]e-2'?(^("°)-*^'"l"»l) g-^^t^^")! (47) 
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where 

SkH^)] = Ygj rfx{(V^(x))^ + 2A^ (ixpV^-W - ^)} (48) 
is the action of 2D Liouville Field Theory (LFT). 



The careful treatment of the quantum LFT in the case = (see for review |T^) enables 
one to find the conformal weights A (e"^'^''') = q{Q — q)-, where Q{g) is the "background 
charge", obtained by imposing conformal invariance of 5*0 [v?]- The authors of work have 
related the value A (e"^'^'^) to the critical exponent Tan{q) in the scaling dependence of the 
average inverse participation ratio (^4|) 

(Q(g,i?)) (49) 



Despite the multiscaling exponent Tan(q') has been computed in the general framework 
of Conformal Field Theory (CFT) few years ago, from our point of view, the geometrical 
interpretation of the multifractal behavior in the model has not yet been cleared up. A 
more physical approach put forward in |0 exploits an analogy between this model and the 
problem of directed polymers on a Cayley tree. This analogy is supported by the fact that 
in both cases the correlation functions grow logarithmically with the distance. For directed 
polymers it is the correlation function of the random potential defined on the tree vertices 
that scales logarithmically with the ultrametric distance (i.e. distance along the tree). The 
same logarithmic behaviour occurs in 2D Gaussian Field Theory ^j. 

We adopt a different point of view. Let us notice that the tree structure (conjectured 
by CChamon et al) emerges quite naturally from the Liouville field theory treated at a 
semi-classical level. Our idea is as follows. Indeed there is not just one saddle point solution 
but a whole orbit of solutions parametrized by SL{2, IR). If one further assumes that the 
integration has to be performed not over the whole group but only over a subgroup (for 
instance Tg), one recovers quite naturally the model defined in section |I|. 

Our starting point is the semi-classical {g 0) treatment of ( ^Tf ) (a similar approach 
can be found in |ll9[|). Using a saddle point method, one is led to the classical equation 

Aif - 2qg6\^ - xo) + 2Xg\^\^'' e'^^ = (50) 
which gives, after integration 

X = q- F/2g (51) 

where 




is the magnetic flux. We then introduce the shifted field 

^(x) = ip{x) - In Ixl'^ (52) 

and taking into account that 
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A(^(x) = Aip{x) - 2nk5^{x) 



(53) 



we end up with the following equation 

+ 27rk5^{x) - 2q g5^{yi - xq) + (2g c/ - F)e~'^'^ = (54) 



The most general solution of (|5^ (away from singularities) in Euclidean space of complex 
coordinate z can be written as follows [|l5l 



e-'- = ^ d.A{z)ckB{z) ^ ^^^^ 



\F-2qg\ {i ^ ^A{z)B( 



z. 



where A{z) and B{z) are correspondingly holomorphic and anti-holomorphic functions of z 
and e = sign(F — 2qg). The semi-classical treatment assumes g to be small, hence in order 
to have real if in Eq.(|55D we should put e = 1, i.e. q <^ F/2g. The relevant solution of (pSj), 
compatible with the singularities, reads 

= 4(A: + 1)^ {z-zf 

F {l + {zzY+^f ^ ' 



The normalization condition of the wave function is then satisfied only if F = 47r(A; + 1). We 
would like to stress that for = we here recover the critical value 47r of the magnetic flux: 
uniqueness of the ground state wave function holds only below this value. It also should be 
mentioned that our analysis does not depend on a peculiar basis of eigenf unctions and the 
results presented here can be extended to wave functions of the form 'il){z) = Pk{z)e~'-^^^\ Pk 
being a polynomial of degree k. It is noteworthy that (p6D is an algebraically decaying wave 
function, what is, following p|, a signature of the existence of prelocalized states. 

Using the fact that the Liouville field is not exactly a scalar but varies under holomorphic 
coordinate transformations z ^ w{z) as 

^{z)^^{w{z))-\n\w\z)l (57) 



one can check that the set of solutions (|56|) is invariant under the following family of trans- 
formations, parametrized by the group PSL(2,C): 



z ^ wu{z- a, 6, c, d) = I (58) 



The orbit of (pci{z) is then given by 



ipci {z; a, b, c, d) = 
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Up to redefinition of the measure dr on PSL{2,€), we restrict the domain of integration to 
P^L(2,IR). 

Due to the angular symmetry of (0), we take points of the form z = i~ p (p E IR), and 



following [|T^ we rewrite (^4])- 



47r(fc + l)^ 



"Det 



Let us denote 



a b 
c d 



and Un 





-1/2 


/ ^ 






JPSL{2,R) 




fl)/2 









p-{k+l)/2 



-2qtpcl{p,T) 



(60) 



(61) 



then we can rewrite (^) as follows 



/•S'fcM JpSLi2,R) 



(62) 



where we have got rid of irrelevant factors and the function J(p, q) reads 



-^(P) Q') = / [2 cosh(i(i, rz/p(i))] 

JP5L(2,R) 



(63) 



Instead of summing over the whole group PSL{2, IR), we restrict the sum over a discrete 
subgroup, Tg in our case. Even if this discretization of the saddle manifold has no evident 
physical justification, we believe that the model obtained yields interesting results. It leads 
to consider the so-called Poincare series (see for review") H, defined as follows 

(64) 



I{p,(l) = H{i,Up{i),q) = [2cos\id{i,TUp{i))\ 



As shown in ^6], the series H does not converge for q < q* with q* depending on Tq only 
(the analysis of the previous sections show that we roughly may set g* ^ go/2). A new 
dependence on p occurs only if H does not converge, we will therefore consider only this 
regime. We must introduce in this case a cut-off N to regularize the series, and finally study 
asymptotics of the finite sum 



\2 cosh. d{i, TV p{i))] 



(65) 



This Poincare series has the same asymptotic properties as the one that defines our model. 
In particular it will exhibit a multifractal behaviour in M. However what really matters for 
a physical system is the multifractal behaviour under transformations parametrized by p. 
We therefore have to relate the behaviour in the group manifold to the behaviour in the real 
space. This will be achieved through a renormalization transformation of the form 
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%(p, = C'p''V(Ar,p)(l, q) 



(66) 



Appendix provides a heuristic derivation which gives 

lM{p,q) = Cp-^'+'^'^'/^Ij^_^o^, {l,q) (67) 

Comparing (p7|) and (|6lD gives /^/'(l,?) ~ Z{2q,Af). Therefore 

/_^(p,g) = Cp-('=+i)i°2/72(2g,Ar+ ^Inp) (68) 

7 

This relation allows to extract the scale dependence in p for a given cut-off A/". Using the 
asymptotics of Z obtained in previous sections yields 

^2qk-2q^ip)\ ^ 2g-(fc+l) ln2/7 -^t(2<7) 



with the notations of previous sections. After integrating over the whole domain we arrive 
at the final expression for the multifractal exponent r^niq)'- 

\n(Q(q,R)) , J (A; + 1)^2 \ , 

rUq) = - lim = 2(g - 1) 1 - ^ q] for g < go/2 70 

-R^oo mil y 7 / 

with go defined in (^). The regular term 2(g — 1) corresponds to the one obtained in |^ 
for (yf ^ 0. Multifractality of the wave function is induced by the quadratic term, which 
is directly related to geometric properties of the saddle point hyperbolic manifold (target 
space), and holds in absence of any random potential in this target space. The regime 
g > go/2 is not affected by these geometric properties. 



IV. CONCLUSION 



The wave function ipk{^) introduced in (^) belongs to the general class of exponential 
functionals of free fields. Such functionals appear in several physical contexts. 

1. The square of the wave function (for A; = 0) may be interpreted as the equilibrium 
Gibbs measure in the random potential fix). In the ID case the problem was first studied 
A rather deep and complete analysis of this problem was recently presented in p5 



m 



2. Exponential functionals of free fields play an important role in the context of one 
dimensional classical diffusion in a random environment. Their probability distribution 
controls the anomalous diffusive behaviour of particles at large time . 
the study of disordered samples of finite length [^p9 
discussed in [051. 



They also arise in 
Some mathematical properties are 



3. In the context of one dimensional localization in a random potential such functionals 
arise in the study of the Wigner time delay [^]. 
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4. The function ipki^) is the ground state wave function of 2D Dirac fermions in a 
random magnetic field with B = A(p. The multifractal behaviour first conjectured in has 
been recently confirmed by an independent investigation based on renormalization group 
method |2^. The scenario of multifractality which is presented here relies mainly on a 



geometric approach to a semiclassical quantization scheme of the Liouville field theory. The 
fact that tree like structure emerges quite naturally in our consideration is an interesting 
feature which obviously deserves further investigation. The multifractality in our approach 
appears as a by-product of isometric embedding of a Cayley tree in the hyperbolic plane. 
The objects which possess multifractal behavior are the moments of the partition function 
defined as sums over all vertices of a Cayley tree isometrically embedded in the hyperbolic 
plane where each vertex carries a Botzmanm weight depending on the hyperbolic distance 
from the root point. No randomness is imposed in the model. 

From the mathematical side our work reveals some interesting links between the theory 
of automorphic functions, invariant measures and spectral theory. We hope to return to 
these problems in a forthcoming publication. 
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APPENDIX A: 

Let us prove that the function u{w) = (1 — wwY f{w) where f{w) is defined in (|36|) has 
the following properties: 

• At all centers w = Wc oi zero-angled triangles tesselating the Poincare hyperbolic unit 
disc u{wc) = 1; 

• The function u{w) has local maxima at the points Wc- 

I. The proof of the first statement implies the proof of the fact that the function u{w) is 
invariant with respect to the conformal transform w^-^\w) of the unit Poincare disc to itself 
where 

. ill^, (Al) 

WWo — 1 

and wq is the coordinate of any center of zero-angled triangle in the hyperbolic Poincare 
disc obtained by successive transformations from the initial one. 

Hence, it is neccessary and sufficient to show that the values u{w = 0), u(w = 
u (w = |e*'^/^^ and u (w = |e*^'^/^^ coincide. Then, performing the conformal transform and 
taking wo = | — f , ^e^'^^^, |e*^'^/^|, we move the centers of the first generation of zero-angled 
triangles to the center of the disc w^^\ Now we can repeat recursively the contruction, i.e. 
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find the new coordinates of tlie centers of the second generation of zero-angled triangles in 
the disc w^^^ and compute the function u{w) at these points, then we perform the conformal 
transform w^'^Hw^^'^) and so on... 



We have at the point w = 0: 



8/3 



while at the point w 



■| the function u{w) can be written in the form 



u w 



1 - ie— 



i7r[l/2+j/{2V3)h 



8/3 / 



(A2) 



Let us use the properties of Jacobi 6'-functions: 



Taking into account 



u w 



we can rewrite 
2^c 



in the form 



34 



(^33/4^)8/3 Ql^^Q^ gi7r[l/2+i/(2v^)]3^ 



keN 
XeR 



8/3 /3 



^'1(0, e 



«7r(l/2+iV3/2)+«7r 



8/3 



(A3) 



(A4) 



w = = 1. 



Thus, u{w = 0) = u (i 

At the point u (^w = ^e^^^^^ we have 



u w = -e 
2 



i7r/6 



2^0 



8/3 



8/3 



1 - 



In the same way we can transform the function u (w = |e*^'^ 



/6 



u\w = -e^^"/^ 



2^0 



^;(0,e^^[-i/2+^v^/2]) 



8/3 



1 - - 

4 



^'1(0, e*'"[l/2+jV3/2]-j7r-) 



8/3 



(A5) 



(A6) 



The transforms (|A4|) -(|A6|) complete the proof of the part I. 

II. Let us prove that the function u{w) has local maxima at all centers Wc of zero-angled 
triangles tesselating the Poicare hyperbolic disc. Actually, the function f{w) by construction 
gives a metric of some discrete subgroup of the group of motions of Poicare hyperbolic disc. 
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Hence the function f{w) cannot grow faster than the isortopic hyperbohc metric (1 — wwY 
and the following inequality is valid 

< u{w) < 1 

for all points w inside the unit disc. But we have shown that u{w) = 1 a.t w = Wc what 
means that the function u{w) reaches its local maxima at the points Wc and at all these 
maximal points the function u{w) has one and the same value u{wc) = 1. The part II is 
proved. 

APPENDIX B: 

Our aim is to extract explicitly the p-dependence of the truncated series (|65D and to 
connect it to J(l, g). More precisely we are looking for a renormalization transformation of 
the form 

IAp. q) = Cp'^Ix'ixA^, q) (Bl) 

Using the correspondance (up to the volume of S0{2)) between PSL{2, IR) and 7Y^, we 
interpret the shift r — > TUp as a change of hyperbolic coordinates — see fig 0. 




FIG. 7. Change of coordinates in hyperbolic Poincare disc 

Note that the expression (|65|) does not depend on the particular representation of the 
hyperbolic 2-space, since the hyperbolic distance is invariant. The only one requirement is 
to define a compatible action of Tq in the space under consideration. We use for conveniency 
the unit disc representation whose center is the image of the point i'p{i) where Vp is defined 
by Eq.(^) in the Ti^ representation. For shortness the generic element of Tg is labelled by r 
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independent on the representation. We parametrize the point TVp{i) by its hyperbohc polar 
coordinates {d\, a]) with the origin at the point i and by "shifted" hyperbohc coordinates 
(c/2, O-D with the origin at the point I'pii). Note that dp = d{i, Vpii)) ~ (A; + 1) Inp, and the 
following "triangle equation" in hyperbolic 2-space holds: 

cosh (i[ = cosh cosh dp + sinh dl^ sinh dp cos (B2) 

In order to extract the scaling conjectured in (^), we make an approximation which consists 
in neglecting fluctuations of Wn{d). In this approximation we can sum over the generations 
n and the angles aj^ of the vertices within each generation (1 < jn ^ 3 x 2"^^). Namely 
{d2,a.2) = ipfn^aj^). Thus one has 

M 3x2"-i 

Ij^f{p,q) = '^[2cos]i'yncos]idp]~'^'' ^ (1 + tanh7ntanh(ipCOSQ;j„)~^^ (B3) 

n=l j„=l 

Assuming that aj^ are uniformly distributed, we get for n ^ 1 the following expression 



3x2 



n-l 



3 X 2"-i /-s^ da 



^ (l + tanh7ntanhdpCosajJ 2'?^ — / ——. (B4) 

Itx Jo (1 + tanh7ntanh dp cos a j^^ 

which leads to the asymptotic behavior: 

3x2"-i r const for g < 1/4 



lim 2 " (1 + tanh7ntanh(ioCOSQ;,„) 

n— >cxD ^ — ' ^ •' 

jn = l 



>oo 



(B5) 

g2dp{2g-5) g > 1/4 



As justified hereafter we consider as relevant only the case q < 1/4. Using that for n ^ 1 
and p ^ 1 one has 2 cosh 7/1 cosh dp ~ cosh(7n + dp) we can rewrite Eq.( P3|) as follows 

Ij^{p, q) = Cp-^^^^)'^^h ^ 2("+'^-/^)[cosh(7n + dp)]-^'^ (B6) 

71=1 

Performing the shift n = n + dp/ 7 we get finally 
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This expression fulfills the condition (pl|) . We assume that this renormalization also holds 
for the full function /^(p, q)- 

Let us pay attention to some contradiction between (^) and (|B5| ) raised by the set of 
successive approximations of (|63D which however is irrelevant for our final conclusions about 
multifractality. The equation ( |55D shows that if the integral over PSL{2, IR) converges, it 
should not depend on p. Using the Poincare series, the convergence of (|63|) occurs for q > q*. 



For g > 1/4 and q > q* the p-dependence shown in (|B5|) should then cancel by summing 
over all n. The discrepancy between ( |63D and (p5| ) appears for q G [1/4, g*]. First of all 
we should note that the interval [j,q*] is numerically small (following previous sections we 
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have g* ~ go/2 — 0.4) and is nonuniversal, i.e. depends on the particular choose of the 
subgroup under consideration. Moreover, both the threshold g = 1/4 and the asymptotics 
( [B5| ) depend on the distribution of aj^ and we believe that more careful treatment of angle 
dependence in (|63D would allow disregard the region [1/4, g*]. 
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